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■ a scheme using a dipole force barrier to select the coldest àtoms from a magnetically trapped atom 

cloud. The dipole and magnètic potentials create a local minimum which traps the coldest àtoms. 
A unique advantage of this technique is the sharp cut-off in the velocity distribution of the sample 
of selected àtoms. Such a non-thermal distribution should prové useful for a variety of experiments, 
including proposed studies of atòmic tunneling and scattering from quantum potentials. We show 
that when the rms size of the atom cloud is smaller than the local minimum in which the selected 
àtoms are trapped, the velocity selection technique can be more eíRcient in 1-D than some common 
techniques such as evaporative cooling. For example, one simulation shows nearly 6% of the àtoms 
retained at a temperature 100 times lower than the starting condition. 



We discuss a velocity selection technique for obtaining cold àtoms, in which all àtoms below a 
certain energy are spatially selected from the surrounding atom cloud. Velocity selection can in 
some cases be more eííicient than other cooling techniques for the preparation of ultracold atom 
clouds in one dimension. With quantum mechanical and classical simulations and theory we present 
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I. INTRODUCTION 



Advances in làser cooling of neutral àtoms during the last two dècades have led to ob- 



servations of many new i 
atom-optics technologiesj 



nenomena 




as well as the development of a variety of 



To study some of these effects, such 



as quantum tunneling, and studies in quantum chaos, only require that the àtoms being 



studied be cold in one dimension 



es in quantum 



For such studies, cooling àtoms in only 
one dimension should serve the purpose and can prové more efïicient than cooling in all 
three dimensions. Cooling in one dimension has been previously studied in processes such 
as 'delta kick' cooling, which compresses the velocity distribution of an atom cloud with 
magnètic íield gradients, paid for with an increase in the spatial distribution, thus conserv- 
ing phase space densitv |l'iT[. As well, àtoms have been cooled to the ground state of a one 
dimensional optical lattice ïsl lïfli. 

In this paper, we analyze an efïicient and straightforward method to achieve low one 
dimensional temperatures in neutral àtoms through velocity selection. The technique takes 
advantage of the fact that in an atom cloud with a Gaussian velocity distribution there 
is a large population near zero velocity along any individual direction. By utilizing the 
dipole force to manipulate the motion and spatial position of àtoms, a process dependent on 
the àtoms' own energies, cold àtoms can be separated from their high energy counterparts. 
Manipulation of àtoms with the dipole force is a commou experimental technique and has 
been used in several ways such as atom trappingj2^ reflection of àtoms [^|, quantum 
state engineeringj^l, and a similar selection of low energy àtoms in an optical lattice [3, 2^ - 
We show here that this version of velocity selection can not only provide cold samples of 
àtoms at high efiiciencies, but even pure 1-D quantum states at high efïiciencies when in the 
ultracold regime of ~ WnK. Experimental studies on this technique were performed in this 
laboratory in parallel with this theoretical study j^l- 



II. VELOCITY SELECTION OVERVIEW 



The proposed velocity selection process is the following: àtoms are trapped in a weak 
quadrupole magnètic trap, while a blue-detuned làser sheet sweeps through the trap. The 
repulsive potential due to the blue-detuned làser pushes the colder àtoms (those with less 
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kinetic energy than the làser potential) away from the more energètic àtoms which can 
classically overcome this potential. The result is a cold sample of àtoms spatially separated 
from the original atom cloud. Away from the center of the magnètic trap the magnètic 
potential is approximately separable allowing us to express the combined magnètic and 
optical potential in one dimension. (This approximation could be made exact by using a 
harmònic rather than quadrupole potential.) A schematic of this is shown in íigure 1. The 
combination of the magnètic and dipole potentials results in a potential of 

U = Um + Ud = fJ'B' \x\ + Uq exp 

where /x is the atòmic magnètic moment, B' is the spatial gradient of the magnètic trap, xd 
is the position of the dipole barrier, Wq is the rms width of the dipole barrier, and Uq is the 
peak of the dipole force potential. To make the above potential conservative, the detuning of 
the làser from the atom's resonant frequency must be large enough so that photon absorption 
by the àtoms is negligible. In typical dipole experiments in which conservative potentials 
are required, a detuning oí 6 > lOOOF is suflficient, where T is the natural linewidth of the 
atom. This combination of the magnètic and dipole potentials is shown in íigure 2. Note 
from íigure 2b that the well depth is not equal to the height of the dipole potential barrier, 
f/o, but rather f/e//, which is somewhat smaller due to the íinite width of the barrier. 

After the selection process shown in íigure 1, the selected àtoms íind themselves in the 
potential well bounded on one side by the magnètic trap and on the other by the dipole force 
beam. The atòmic velocity distribution of this selected cloud is a truncated Gaussian, with 
no àtoms having energy greater than the depth of the potential well. It's important to note 
that when an atom cloud is selected at a temperature we define as Tg// = 2Ueff/kB, where 
ks is Boltzmann's constant, this does not describe the rms kinetic energy, but the absolute 
maximum kinetic energy of the hottest àtoms in the sample. Such a truncated distribution 
is required, in particular, for a planned experiment studying the tunneling of ^^Rb through 
a relatively macroscopic barrierQ. It must be known for certain that no àtoms in the 
sample have energy greater than the barrier height, or else they could traverse the barrier 
classically. In a thermal distribution, there would always be the possibility of classical 
traversal by àtoms in the high energy tails of the distribution. Similar considerations apply 
to other experiments, such as proposed studies of quantum potential scatteringj^!- 



[X 
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III. CLASSICAL THEORY AND SIMULATION 



Until the temperature of the atom cloud is in the ultra low regime where quantum effects 
begin to dominate, the velocity selection process can be largely described by classical physics. 
By considering the magnètic and dipole potentials as classical, conservative potentials, the 
dependence of the efficiency of the selection process on various parameters can be studied 
when still in this regime. 

We begin with a description of the 'experiment' and analytic estimates of selection ef- 
ficiency before presenting the results of the classical simulations. We start with an atom 
cloud with Gaussian spatial and velocity distributions of rms radius Tq and rms velocity l'o 
respectively. The initial temperature is Tq = mvQ/kB- 

We model a dipole force beam moving slowly through the atom cloud. When the dipole 
beam is moved slowly enough, whether or not an individual atom is swept up by the barrier 
depends only on whether its total energy is greater than that of the potential created by 
the barrier. Of principal interest is how the efficiency of the selection process depends on 
barrier height for different atom cloud sizes, temperatures, and magnètic field gradients. 
Three separate regimes occur, in which the dependence of efficiency on barrier height differs 
drastically. Which regime the experiment is in depends on the ratio of the initial kinetic 
energy {{KE). = ^mv^) to the variation of magnètic potential energy across the cloud 



The íirst regime {{KE)^ ^ {PE)^) occurs when the rms radius of the atom cloud is 
much smaller than the size of the potential well formed by the magnètic gradient and dipole 
barrier (ro <^ r^). Since the magnètic trap is linear in position, ~ Ueff/{nB'). When 
1^0 very few àtoms are given significant potential energy by the magnètic gradient, 

as most àtoms are always near the minimum of the potential well. Therefore the velocity 
selection process selects àtoms by their kinetic energy only. The requirement for selection 
of an atom in this regime is merely ^mv'^ < Ueff, ov v < Vc = ·\/2[/e///m = ^/kBT^ffJrn. 
We define the efficiency of the process to be the ratio of number of àtoms selected during 
the process to the total number of àtoms in the cloud. In this regime this efficiency is given 
by an error function of the original Gaussian velocity distribution 



((PE),= ^/x5Vo). 




(2) 
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At low barrier heights, that is, for Tg// < Tq, this can be approximated by 



2 (Vr\ l2Teff 



This is the desired regime for the experiment, as the efficiency of the selection process may 
be high even for large temperature ratios. For example, àtoms as cold as lOnK (similar 
to temperatures required for future experiments) could be obtained from an atom cloud at 
IfiK with an efficiency of ~ 8%. This is to be contrasted with typical efficiencies of < 1% 



23, 



for evaporative cooling by a factor of 100 in temperature 

The other extreme occurs when tq ^ {{KE). < {PE)-), that is, when the atom cloud 
size is much greater than the size of the potential well. In this case, even though the initial 
atom temperature, Tq, may be much less than the dipole potential height, a significant 
amount of potential energy is given to many of the àtoms by the magnètic gradient. The 
requirement for selection of an atom now involves both the kinetic and potential energy 
of the atom, Ua + KEa < U^jf. During the selection process, àtoms are being selected 
by their kinetic energy as well as their potential energy, given by their position in the 
potential well. Consider the requirement for an atom of zero kinetic energy to be selected. 
The atom must be close enough to the well minimum that its energy is less than the well 
depth; hbB'x < Ueff, or x < Xc — Ueff/{fiBB') = kBTeff/{2jj,BB'). For an atom cloud 
with a gaussian spatial distribution then, the efficiency of selection, in terms of the spatial 
component only, is another error function integral: 

1 



riPE 





X^ 


1 exp 


. ^' 0. 



dx (4) 



'2'KrQ 

Again, for low barrier heights such that Tgjj < Tq, this can be approximated as r]pE 



y'2/7r(xc/ro) — ^j2pRTf,f fkB j {2ri^\XBB'\ In conjunction with the square root dependence 
of efficiency on barrier height when only the kinetic energy is being considered, the expected 

3 /2 1 /2 

efficiency, in the low barrier limit, becomes r] ~ P^^fj, where P = kB/iT^ roTTfiBB'). A 
dependence of this sort is less desirable than that of equation (3), as the efficiency falls off 
faster than the final temperature. (Note that in evaporative cooling, efficiency is typically 
roughly linear in temperature.) 

In between these two regimes, the efficiency should be approximately linear with barrier 
height, 7] oc Teff. The exact expression for the efficiency in all cases for this description of 
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the selection process is 

1 

T) — / 

27rvoro J_^, 

Note that since, in the general case, the critical position for selection is Xc = {Ueff — 
^mv"^) / {^bB') , this equation is quite convoluted. This expression has been numerically 
evaluated for various regimes and gives the same qualitative dependences as the actual 
classical simulation. 

Figure 3 shows the results of the simulation for selection efïiciency versus barrier height 
for a very small atom cloud (ro = 25/im) of temperature Tq = bQjiK in a well made of 
a magnètic gradient of B' = 0.5G/cm , and a dipole barrier with an rms radius of 
For comparison, the sòlid line presents the results of equation (3). For these parameters 
(KEi) — 373 (PEi). The graph definitely shows an 77 ~ ^/T^JfJ^o dependence of efïiciency 
on barrier height. As many as 5% of the àtoms may be selected at temperatures as low as 
Te// = òOOnK = To/lOO. 

Next, shown in figure 4, is the efïiciency versus barrier height for a much larger magnètic 
gradient of lOOG/cm, with the other parameters unchanged. Now (KEi) = 2 (PEi). This 
shows much lower efïiciency because of the spatial dependence now present. The accompany- 
ing sòlid line shows the l3T^0 dependence suitable for this regime. Note from the condition 
for the large cloud size limit, ro ^ C/e///(A*B-B'), that a larger initial atom cloud size, as 
opposed to a larger magnètic field, would also create an efïiciency of this dependence. 



2vl 



J Xn 



exp 



X 

2^ 



dxdv. 
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IV. QUANTUM THEORY AND SIMULATIONS 



It is at ultra low temperatures that quantum effects start to become evident. Performing 
simulations using quantum theory allows us to explore these low energy effects not predicted 
by the classical description. Instead of considering a classical Gaussian cloud of àtoms, we 
now start with an initial wavefunction í = 0) given by 



i,{x,t = 0) = ^v„feí = o)p„'/V 



E 

n 



27rxl 



1/4 



exp 



( X OC f 



4x1 



Pl'\Xn)ét'r^ 



where P^i 



1/2 



(6) 



X 

'2r. 



J 



is the weighting due to the overall gaussian distribu- 



tion of the ensemble, ipn{x,t = 0) = 



2lTX 



.1/4 


(x-Xn)'^ 


j exp 


[ 4^0 J 



, and (f)n is a random phase. 



Random phases were given to each component to ensure that there is no phase coherence be- 
tween individual components, making the initial wavefunction a collection of highly localized 
energètic àtoms. To make sure this is the case we analyze the initial momentum distribution 
with a Fourier transform and check that it resembles a fairly smooth gaussian distribution 
of rms momentum width po = h/{2xo). The result is one representation of a thermal 
cloud of àtoms of rms radius ro, and initial kinetic energy {KE)^ = pl/{2m). For these 
simulations we use the mass of ^^Rh. The wavefunction ip{x,t) obeys the time dependent 
Schrodinger equation ih-^'ip{x,t) = {Hq + U{x,t)) il>{x,t) where Hq = — fi?/{2m)d'^/dx'^ and 
U{x,t) is the combined magnètic and dipole potential. The solution to the time dependent 
Schrodinger equation is found by using the split-operator method and evolving the initial 
wavefunction in time by 

»í+Aí 



ip{x,t + At) = exp 



{Ho + U{x,t))dt 



Íj{x,t) 



exp 



—At— 

4m dx'^ 



exp 



--AtU{x,t) 

1 1 



exp 



4m dx'^ 



Í){x,t). 



(7) 



This equation is evaluated in momentum space. A Fourier transform then yields the final 
real space wavefunction ip{x,tf). The final result of the velocity selection is then found 
from \ílj{x,tf)f to obtain the probability distribution. Because of the random phases be- 
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tween each ipniXjtAe 



[xA 



^„Pri(a;n) \'^n{,x,tf)\^ ^ makiug clear the connection of this approach to a more computa- 
tionally intense density-matrix approach. 

Shown in figure 5 is a plot of the total probability distribution after selection 

in a magnètic trap of IQnK/ fim gradient along with a dipole barrier with a 97nK depth and 
5/im rms width which, together, make a potential well of A.GnK. Also shown is the potential, 
U{x,tf), showing a local minimum at 27 fim in which selected àtoms are trapped. The 
original atom cloud has an rms radius of tq = 8/im, while the rms width of each individual 
wavepacket, 'ipn{x,t = 0), is Xq = 0.07/im, corresponding to a temperature of Tq = 292nK 
(a reasonable starting point after techniques such as delta-kick cooling)jï^. The result of 
the selection is a thermal atom cloud which has expanded in the weak magnètic field due 
to its kinetic energy. At the local minimum is the selected portion of the wavefunction, 
representing 5.8% of the original wavefunction and having a width of ~ 2jj,m. The energy 
of the selected atom cloud, ips{x), is found by projecting out ÍJs{x) and calculating the 



can also be written as 



expectation vàlues of its potential and kinetic energy. Potential energy is given by 



where the low point of the well Umin is subtracted so that the result is with respect to the 
local minimum. For the kinetic energy, a transform into momentum space, ips{p), is made 
and the expression 



is evaluated. The results for the wavepacket shown in figure 5 are a potential energy of 
l.GnK and a kinetic energy of l.lònK. Close to the minimum of the potential well, the 
shape is very close to a harmònic potential. If a harmònic approximation is made, one 
obtains a trap frequency (for ^^Rb) of 280Hz corresponding to a ground state energy of 
2.1nK, and first excited state energy of. S.lònK While the selected wavefunction shown 
in figure 5 certainly looks like the ground state of a harmònic oscillator, the diíference in 
potential and kinetic energy in the cloud and its diíference in total energy from the expected 
ground state energy points to some slight anharmonicity in the potential well as well as a 
slight mixture with the first excited state. 

Shown in figure 6 is a similar plot with the same atomic-cloud parameters as for figure 
5, but here a larger dipole potential of 107nK creates a deeper well of C/g// = 9.8nK. The 
result is a slightly larger selected wavefunction no longer resembling the ground state of the 
well. 

The quantum simulation also yields a prediction for the expected selection efficiencies. 
By integrating over the probability distribution for the selected cloud, \'iljs{x)f , the efficiency 
of selection for the A.6nK and 9.8nK well depths shown in figures 5 and 6 were found to 
be 5.8% and 8.5% respectively. By varying the height of the dipole barrier, a graph of the 
efficiency versus well depth may be constructed with the quantum theory for comparison 
to the classical predictions. Figure 7 shows the dependence of the efficiency on the well 
depth, Ueff for the same parameters as used for figures 5 and 6 but with a varying potential. 
Also shown, in dashed lines, is the efficiency versus well depth as calculated by the classical 
simulation for the same parameters. In the inset are shown the results of the quantum 
simulation for very small well depths. The quantum simulation shows general agreement 
with the classical simulations. For low dipole barriers however, the quantum simulation 




(8) 




(9) 
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produces an efïiciency curve with some step structure, indicating energy gaps for which an 
increase in barrier height does not show a corresponding efRciency increase. 

V. CONCLUSION 

We have proposed a method of not only achieving ultra-low one dimensional temperatures 
in samples of àtoms, but also of creating speciíic quantum states and nonthermal mixtures 
in a potential well. For very small clouds of < 50/im a large number (~ 5%) of àtoms can 
be selected with an efRciency which varies as the square root of the final temperature. This 
is much more economical than techniques whose efíiciency drops linearly with temperature. 
The technique also guarantees a maximum energy set by the parameters of the experiment. 
This sharp cut-off in the velocity distribution is useful when performing experiments in which 
the energy of each atom must be below a certain limit, where a Boltzmann distribution 
would be unsatisfactory. Our theory and simulations show that there are experimentally 
accessible parameter regimes where velocity selection can be more efíicient in terms of final 
atom number than existing cooling methods at similar temperatures. Àtoms at temperatures 
of < 20nK are shown to be in non-classical states, but as higher temperatures are selected, 
the resultant cloud is increasingly classical. These characteristics of the velocitv selection 
technique provide the motivation for present experimental work on the method |25l|. 
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Figure 1: Schematic of velocity selection process. Àtoms are first placed in a shallow 
magnètic trap. Moving dipole potential then sweeps cold àtoms up the magnètic potential 
away from more energètic àtoms. 

Figure 2: Shape of combined potential of quadrupole magnètic fields and blue detuned 
dipole force làser beam. Close-up shows effect of gradient of magnètic potential and width 
of dipole beam on the well delpth. 

Figure 3: Simulated efíiciency versus barrier height for a 25//m, òO/iK atom cloud in a 
O.SG/cm magnètic potential. For low C/g// (<S kBTç,), the efíiciency can be approximated 
by the analytic expression rj = ^y 2Tef f / (ttTq) , shown by the sòlid line. This favorable 
dependence is seen for parameters such that fisB'ro <^ Ueff, illustrating the regime where 
negligible potential energy is given to the atom cloud by the magnètic trap potential. 

Figure 4: Simulated efíiciency dependence graph for a 25//m, òO/iK atom cloud in a 
lOOG/cm magnètic potential. With these parameters, ^bB'tq ^ C/g// and efíiciency now 
has a dependence on both the velocity and spatial components of the atom cloud, resulting 
in a dependence at low well depths which can be approximated by 77 = PTeff^^"^ which is 
shown as the sòlid line for comparison. 

Figure 5. Probability distribution, |\í'(a;, íj)|^, of a thermal atòmic cloud after velocity 
selection. The selected cloud, ipsi^), is towards the left of the plot in a 4.6nK local minimum, 
spatially separated from the rest of the more energètic àtoms remaining in the center of the 
trap. 5.8% of the original atom cloud was selected to form this wavepacket having potential 
and kinetic energies of 1.6nK and l.lònK respectively. The initial atom cloud had a rms 
size of and an rms temperature of 292nK making a ratio of final to initial temperature 
Tf/To ~ 1/100. Note the resemblence of the selected probability distribution to that of a 
single pure quantum state in this harmonic-like potential. Also shown is the combination of 
the magnètic and dipole potentials showing a shallow well where the selected wavefunction 
is resting. 

Figure 6: Plot of for a deeper potential well of 9.8nK. Àtoms were selected 

into the local minimum here with an efíiciency of 8.5% and a ratio of final to initial tem- 
perature oíTf/Tç) ~ 1/30. At this deeper selection potential more than one bound state is 
observed as shown by the structure in the selected wavepacket. 

Figure 7: A plot of efficiency of selection versus the potential depth, Ueff calculated 
from the quantum simulation (sòlid line) as well as from the classical simulation (dashed 
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line) for comparison. Insets show the selected wavefunctions, ipsi^), from figures 5 and 6 
and their relative positions on the efïiciency curve. 



13 



